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We study roughness probability distribution functions (PDFs) of the time signal for a critical 
interface model, which is known to provide a good description of Barkhausen noise in soft ferro- 
magnets. Starting with time "windows" of data collection much larger than the system's internal 
"loading time" (related to demagnetization effects), we show that the initial Gaussian shape of the 
PDF evolves into a double-peaked structure as window width decreases. We advance a plausible 
physical explanation for such structure, which is broadly compatible with the observed numerical 
data. Connections to experiment are suggested. 
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I, INTRODUCTION 

The probability distribution functions (PDFs) of crit- 
ical fluctuations in assorted systems have been the sub- 
ject of much recent interest, mainly stemming from the 
realization that they exhibit a remarkable degree of uni- 
versality diii. 

The roughness w 2 of a fluctuating interface with N 
elements is the position-averaged square width of the in- 
terface height above an arbitrary reference level : 
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w 2 = N- x J2( h i- h Y 



(1) 



where h is the average interface height. The finite-size 
scaling of the first moment of the roughness PDF gives 
the roughness exponent £ 0: 



(w 2 (L))~L 2 < , 



(2) 



where angular brackets stand for averages over the en- 
semble of allowed interface configurations, and L is some 
finite linear dimension characterizing the system in study. 
The width PDF P{wi) for correlate^ systems at critical- 
ity may be put into a scaling form 
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<5>(z) = (w 2 ) P(w 2 ) , Z = W 2 /{w 2 ) , 



(3) 



i.e., the scaling function $(z) is expected to depend only 
on the scaled width w 2 /(w 2 ). In other words, the size 
dependence must appear exclusively through the aver- 
age width (w 2 ). Comparison of experimental or simula- 
tional data to specific analytical forms, whose suitability 
to the description of the case at hand has been antic- 
ipated by physical arguments, usually results in good 
agreement. Thus, the PDF of voltage fluctuations in 
semiconductor films was fitted very well by that of perfect 



Gaussian 1/f noise 0; simulational data for the single- 
step model of deposition-evaporation, by the PDF of a 
random-walk process (the latter corresponds to per- 
fect Gaussian 1/f 2 noise, or Wiener process 0). Further 
progress was made possible via the analytical evaluation 
of roughness PDFs for generalized Gaussian noise with 
independent Fourier modes ( i.e. 1//" noise with gen- 
eral, continuously-varying a) 0| ■ Consideration of the 
scaling properties of height-height correlation functions 
and their Fourier transforms implies that 



a = d + 2( 



(4) 
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where d is the interface dimensionality and £ is defined 
in Eq. © . Eq. Q is valid provided that C > 0, i.e. 
a > 1 in the present case where the "interface" is a time 
series (see below). 

In previous work we applied the ideas outlined 
above, to investigate the interface roughness PDFs of a 
single-interface model which has been used in the descrip- 
tion of Barkhausen "noise' ' (BN) 0GSEHI3. This is 
an intermittent phenomenon which reflects the dynamics 
of domain- wall motion in the central part of the hystere- 
sis cycle in ferromagnetic materials (see Ref. for an 
up-to-date review). By ramping an externally applied 
magnetic field, one causes sudden turnings (avalanches) 
of groups of spins. The consequent changes in magnetic 
flux induce a time-dependent electromotive force V(t) on 
a coil wrapped around the sample. Analysis of V(t), as- 
sisted by suitable theoretical modeling, provides insight 
into both the domain structure itself and its dynamical 
behavior. It has been proposed that BN is an illustration 
of "self-organized criticality" El El El, 

in the sense 

that a broad distribution of scales (i.e. avalanche sizes) 
is found within a wide range of variation of the external 
parameter, namely the applied magnetic field, without 
any fine-tuning. The interface model studied here in- 
corporates a self-regulating mechanism, in the form of a 
demagnetization factor. 

We have shown that the demagnetizing term is 
irrelevant as regards interface roughness distributions, 
with the conclusion that in this respect the behav- 
ior of self-regulated systems is in the same universal- 
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ity class as that of the quenched Edwards- Wilkinson 
model 0, 0, 0, I2D], at criticality (i.e. at the inter- 
face depinning transition). 

However, when one considers the time series of in- 
termittent events which characterizes BN, it is known 
that the demagnetizing term is responsible for the in- 
troduction of short-time negative correlations in the 
model (such correlations are observed in experiments as 
well) |2|. The question then arises of whether a corre- 
sponding signature of self-regulation will be present when 
the roughness distribution of the time sequence of BN 
events is examined. Since the traditional data acquisi- 
tion method in the study of BN is exactly via the time 
series of induced voltages, an investigation along these 
lines may establish useful connections between observa- 
tional data and the basic physical mechanisms underlying 
BN. 



II. MODEL INGREDIENTS AND DYNAMICS 

We use the single-interface model introduced in Ref. @ 
for the description of BN. In line with experimental pro- 
cedure, the external field H acting on the sample is as- 
sumed to increase linearly in time, therefore its value is 
a measure of "time". We consider the adiabatic limit of a 
very slow driving rate, thus avalanches are considered to 
be instantaneous (occurring at a fixed value of the exter- 
nal field). In this simplified version, a plot of V(t) against 
t consists of a series of spikes of varying sizes, placed at 
non-uniform intervals. Generalizations for a finite driv- 
ing rate may be devised fiH I2IL l22l | , but will not concern 
us here. 

Simulations are performed on an L x xL y xoo geometry, 
with the interface motion set along the infinite direction. 
Since we are interested in fluctuations of the Barkhausen 
signal in time, we keep geometric aspects at the simplest 
level, i.e. L y = 1 (system dimensionality d = 2, interface 
dimensionality dl = 1). Periodic boundary conditions are 
imposed at x = 0, L . 

The interface (180-degree domain wall separating spins 
parallel to the external field from those antiparaiiei to it) 
is composed by L discrete elements whose x coordinates 
are Xi — i, i — 1,...,L, and whose (variable) heights 
above an arbitrary reference level are hi. The simulation 
starts with a flat wall: hi — for all i. 

Each element i of the interface experiences a force given 
by: 

fi = u(xi, hi) + k [hi+i + hi-i - 2 hi] + H e , (5) 

where 

H e = H — r)M . (6) 

The first term on the RHS of Eq. (JjJ represents quenched 
disorder, and is drawn from a Gaussian distribution of 
zero mean and width R; the intensity of surface tension 
is set by k, and the effective field H e is the sum of a 



time-varying, spatially uniform, external field H and a 
demagnetizing field which is taken to be proportional to 
M = (1/L) J2i=i hi-, the magnetization (per site) of the 
previously flipped spins for a lattice of transverse width 
L. Here we use R = 5.0, k = 1, rj = 0.005, values 
for which fairly broad distributions of avalanche sizes are 
obtained @, Efl El E3 . 

The dynamics goes as follows. For fixed H, starting 
from zero, the sites are examined sequentially; at those 
for which fi > 0, hi is increased by one unit, with M 
being updated accordingly; the corresponding new value 
of u is drawn. The whole interface is swept as many times 
as necessary, until only sites with fi<0 are left, which 
marks the end of an avalanche. The external field is then 
increased until fi — for at least one site. This is the 
threshold of a new avalanche, which is triggered by the 
update of the site(s) with fi = 0, and so on. 

The effect of the demagnetizing term on the effective 
field H e is that at first it rises linearly with the applied 
field H, and then, upon further increase in H, saturates 
(apart from small fluctuations) at a value rather close 
to the critical external field for the corresponding model 
without demagnetization 0, 0] . 



III. TIME SERIES: CORRELATIONS AND 
ROUGHNESS 

As explained above, owing to the assumed linear in- 
crease of applied field with time (in analogy with exper- 
imental setups), we shall express time in units of H as 
given in Eqs. (JJ and We have generated time series 
of BN, with £>(10 4 - 10 5 ) events. Steady state, i.e., the 
stabilization of H e of Eq. @ against external field H, 
occurs after some 200 events, for the range of parameters 
used here. Though we used only steady-state data, it was 
noted that inclusion of those from the transient does not 
appreciably distort any of the quantities studied. 

In experiment, the integrated signal f At V(t) dt is pro- 
portional to the magnetization change (number of up- 
turned spins) during the interval At. In the adiabatic 
approximation used here, a box-like shape is implicitly 
assumed for each avalanche (i.e. details of the internal 
structure of each peak, as it develops in time, are ignored, 
on acount of its duration being very short), thus the in- 
stantaneous signal intensity (spike height) is proportional 
to the corresponding avalanche size. 

We sample the fluctuations of the signal along succes- 
sive "windows" of equal time duration W, each containing 
many spikes. Each window is divided into equally-spaced 
bins of size 5; the signal intensity associated to each bin 
is the sum of the sizes of all avalanches which occurred 
within that bin. The roughness W2 of the signal on a 
given window starting, say, at t = 0, is given by 

1 W/S 

^ = 1^E(^-^) 2 > £ V(t), (7) 

' *=1 te\(i-X)5,i S] 
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Figure 1: (Color online) Double-logarithmic plot of probabil- 
ity distribution, P(WT), of waiting times (lower curve), and 
accumulated distribution, P aC c(WT) = f Q T P(t) dt (upper 
curve). L = 400. 



where V is an average of V(t) over the whole window 
span W. 

As the signal is intermittent, there are significant pe- 
riods (waiting times, henceforth referred to as WT) of 
no activity at all. Such quiet intervals must be properly 
accounted for in the statistics of fluctuations, hence care 
must be taken when setting up the bin size 6. 

We have examined WT distributions, for varying lat- 
tice widths L = 200, 400, 800. In Fig. Q] (lower curve) 
we display a double-logarithmic plot of the probability 
of occurrence of assorted WTs for L = 400, against WT, 
sampled over 8 x 10 6 events. The distribution is gen- 
erally rather flat, apart from (i) a sharp cutoff at the 
high end (related to the finite cutoff in the avalanche 
size probability distribution, see the discussion of load- 
ing times below), and (ii) a number of peaks concen- 
trated in a somewhat narrow region corresponding to 
10~ 5 < WT < 10~ 4 . The latter are associated to very 
frequent and small, spatially localized (i.e. non-critical) 
events involving typically N = 1 — 10 sites [11]. This 
is easy to see by recalling from Eqs. JU and @ that, 
since the demagnetization term keeps H e approximately 
constant, a small avalanche with N spins overturned de- 
creases the internal field by r]N/L, thus requiring ap- 
proximately the same increase in external field in order 
to bring the system back to criticality. We have checked 
that the peaks move consistently with this argument, i.e. 
their horizontal position is shifted leftward by a factor of 
log 10 2 for each doubling of L. 

Upon consideration of integrated WT distributions 



(upper curve in Fig. QJ, we decided to set the bin size 
S = 10~ 5 (for L = 400, the system size for which most 
of our calculations are done, see below). With such a 
choice, WTs shorter than 5 occur with less than 1% fre- 
quency. This ensures both that inactive periods are not 
wrongly obscured by bursts of activity, and that consecu- 
tive avalanches are rather unlikely to be lumped together. 

At this point, a comment must be made on the con- 
nection of the above results with previous investigations 
of WT distributions in BN. In Ref. it was predicted, 
from a fractal analysis of the ABBM j24[ model which 
describes domain-wall motion via a Langevin equation, 
that P(WT) ~ (WT)~( 2 ~ C ), where c is proportional to 
the external field driving rate. Experimental data in SiFe 
samples are consistent with this [23]. The present case 
of adiabatic driving would then correspond to c — > 0. 
However, it is crucial in the analysis of Ref. |2^ that the 
BN pulse durations be finite (even though they shorten 
accordingly in the c — > limit). Indeed, the result just 
quoted relies on considering the properties of complemen- 
tary sets, both with non-zero fractal dimension (namely 
the time intervals during which there is domain- wall mo- 
tion, versus those of no activity, i.e. WTs). The approxi- 
mation used here, of considering BN pulses as having ex- 
actly zero duration, destroys the connection of our data 
with the conceptual framework in which the power-law 
dependence P(WT) ~ (WT)-( 2 " C ) was found. Though 
in this sense the flat WT distribution found here is most 
likely an artifact of the model, the conclusions extracted 
from the distribution with respect to the choice of 5 re- 
main valid. 

We now turn to the choice of window width W. Re- 
call that real-space properties, e.g. interface roughness, 
of the systems under study benefit from divergence of 
the system's natural length scale, as it self-tunes its be- 
havior to lie close to a second-order (depinning) transi- 
tion @. For such quantities, universality ideas apply, 
so one expects finite-lattice effects to be present only as 
an overall scale factor, e.g. (102) in Eq. © >l Q la u], 
However, in the study of time series for the same sys- 
tems, one must bear in mind that a finite time scale tx 
("loading time") is introduced via the demagnetization 
term Q. This is illustrated in Fig. (similar plots, ex- 
hibiting both simulational and experimental results, can 
be found in Ref. 0) where normalized two-time correla- 
tions (V(t)V(t + T))/(V{t)) 2 - 1 (averaged over t) are 
shown. Therefore, different regimes will be found, de- 
pending on the value of x = W/tl- The limit x ^> 1 is 
expected to reproduce the white noise characteristic of 
uncorrelated fluctuations, for which the roughness distri- 
bution is a pure Gaussian. On the other hand, non-trivial 
effects may arise for x ~ 1. 

Before going further, it must be remarked that tl 
in fact decreases for increasing L. This can be un- 
derstood by recalling that (i) the probability distri- 
bution for avalanche size s g oes roughly as P(s) ~ 
s~ T exp(— s/sq) P. ITnL UlL Il2l|: (ii) the cutoff s scales 
approximately as sq ~ L°- s in the present case of a one- 
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Figure 2: (Color online) Normalized two-time correlations 
(averaged over t) (V(t) V(t + r))/{V{t)) 2 - 1 for system with 
L = 400. Dashed line is fit of data to single-exponential form, 
from which tl = 0.14(1) . 



dimensional interface [lfj|. Thus the maximum waiting 
time tm will vary as tm = V sq/L ~ L~ 0/2 . For L = 400, 
we find so — 4 x 10 4 [12], which explains both the sharp 
drop in the WT distribution at WT ~ 0.5 in Fig. and 
the complete vanishing of correlations at r > 0.5 in FigEl 
In BN studies the connection between lattice-size- 
dependent quantities in simulations, and their exper- 
imental counterparts, becomes especially clear when 
one considers the L-dependent cutoff in the power-law 
avalanche size distribution, and its relation ship to the 
maximum domain size in magnetic samples |lfl |. In the 
present case it should be stressed that finite loading 
times are measured in experiment, under suitable condi- 
tions Thus, we assume that the loading times found 
here are not simply a finite-size artifact of simulations, 
bound to vanish in the thermodynamic limit character- 
istic of real systems. Instead, although we are not in a 
position to propose quantitative comparisons, they must 
correspond to the experimentally-observed ones . 



IV. RESULTS 

By generating many realizations of the roughness w 2 
defined in Eq. <0 for given values of the physical pa- 
rameters, we have obtained the corresponding roughness 
PDFs. The shapes of roughness PDFs found here do not 
usually conform to the generalized Gaussian (1// Q ) dis- 
tributions introduced in Ref. 0, although they display 
certain similarities to the pure Gaussian limit, which cor- 



responds to a = 1/2 in the scheme of Ref. 0. We have 
found it convenient to adhere to conventions used in that 
Reference and related work, namely expressing the PDFs 
in a scaling form, see Eq. Q). 

We first examine the limit x 3> 1. Similarly to l/f a 
PDFs with a < 1 our results in this limit approach a 
S- function shape when expressed in terms of z of Eq. 0- 
The solution, pointed out in Refs. 0, 0, is to use scaling 
by the variance, instead of by the average, i.e. switch to 
the variable 

= W 2 ~ (W 2 ) , . 

The corresponding scaling function will be denoted by 

In Fig. we show results for window width W = 100, 
in terms of y of Eq. JjJJ. While in (a) the demagnetiz- 
ing factor is 77 = 0.005 (thus tl — 0.14 from Fig. [2j, the 
data in (b) correspond to simulations of the same sys- 
tem, with 7] = 0. As explained in Ref. 0, in this case the 
system is kept close to criticality by the following proce- 
dure. We first determined the approximate critical value 
-ffg of the internal field H e of Eq. , by starting a sim- 
ulation with r\ ^ and waiting for H e to stabilize. At 
that point, we set 77 = and repeatedly swept H in the 
interval (7 H£, Hf), 7 < 1, according to the procedure 
delineated in Sec.|n| We have used 7 = 0.9 for the data 
displayed in Fig. El(b). With H° ~ 5.4 for the disorder 
and elasticity parameters used here, data corresponding 
to a window of "width" W = 100 in this case was in fact 
given by the collation of data from ~ W/ (1 — 7) H° = 185 
consecutive field sweeps as just described. Note that, 
within a given field sweep, many non-critical events are 
thus sampled (which would by themselves give rise to a 
non-universal PDF, see below the discussion for narrow 
windows). However, owing to the central limit theorem, 
the result of the collation of many independent segments 
should yield an overall behavior which is essentially Gaus- 
sian. 

One can see that in both cases, a single Gaussian cen- 
tered at y ~ and with variance ~ 1 gives a good fit to 
data, confirming our expectation that demagnetization- 
induced correlations would be essentialy washed away for 
W ^ tl- It is worth mentioning, however, that the un- 
sealed variables tell a slightly different story: for the data 
of Fig. (a) one has (w 2 ) ± a = (127 ± 6) x 10 3 , while 
in (b) (w 2 ) ± a = (6.3 ± 2.2) x 10 3 . Clearly, our data 
would approach a 5- function shape if plotted in terms of 
z defined in Eq. 0. 

Considering now narrower windows, and keeping the 
demagnetizing factor r/ = 0.005, we show data for W = 
10.0, 2.5, and 1.0 in Fig. 0] where we have reverted to 
plotting our results in terms of the variable z defined in 
Eq. 0. This is because it was noticed that, against di- 
minishing x, the scaled roughness PDFs followed a trend 
away from the 5-function shape which was the motiva- 
tion for using the variable y of Eq. © . In order to pro- 
duce an accurate picture of deviations from the Gaussian 
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Figure 3: (Color online) Scaled roughness distributions ^(y) 
of time series, for y of Eq. Jgjl. L = 400, window width 
W = 100. (a): demagnetizing factor r\ = 0.005 (tl — 0.14), 
6x 10 3 samples. Dashed line is Gaussian fit to data with mean 
at y = 0.04(1), width a = 0.96(1). (b): demagnetizing factor 
77 = (see text), 2.1 x 10 4 samples. Dashed line is Gaussian 
fit to data, with mean at y = —0.07(1), width a — 0.98(1). 



Figure 4: (Color online) Scaled roughness distributions $(2) 
of time series, for z of Eq. 0. L = 400. Window width 
W = 10 (triangles, 1.2 x 10 5 samples), 2.5 (squares, 1.2 x 10 5 
samples), and 1.0 (crosses, 5.7 x 10 5 samples). 



limit, we have generated a much larger number of sam- 
ples (O(10 5 )) than for W = 100. 

Before analyzing the shapes exhibited in Fig. 31 it is 
instructive to check how the demagnetization term influ- 
ences the roughness PDFs in the narrow- window limit. 
In Fig. 0the scaled distributions for W = 10 are shown, 
both with and without demagnetization. The shapes of 
PDFs are clearly rather distinct from each other, high- 
lighting the relevance of demagnetization effects in this 
limit. For 77 = the distribution peaks at z ~ 0.15 and 
decays very slowly afterwards. As mentioned above in 
connection with the data of Fig. (b), this reflects the 
non-universal statistics of non-critical events which our 
calculational method for 77 = inevitably includes. The 
difference relative to that case is that for W = 10, each 
roughness sample is the collation of only ~ 19 consecu- 
tive field sweeps. The corresponding results show that, 
in contrast to W — 100, here one is outside the range of 
applicability of the central limit theorem. 

From now on we shall only deal with r\ ^ 0. Even 
though W = 10.0 corresponds to x ~ 70, it is clear from 
Fig. that a secondary peak is evolving, i.e. a signifi- 
cant distinction is emerging with respect to the simple 
Gaussian picture found for larger W, Data for W = 5.0 
(not shown) are virtually identical to those for W — 10.0. 
While a secondary peak still shows up for W = 2.5, data 
for W = 1.0 display only a single maximum (however, 
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Figure 5: (Color online) Scaled roughness distributions $(z) 
of time series, for z of Eq. 10, with and without demagnetiza- 
tion. L = 400, window width W = 10. Triangles, r) = 0.005, 
1.2 x 10 5 samples; circles, r] = 0, 2.1 x 10 5 samples. 
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Figure 6: (Color online) Scaled roughness distribution $(z) 
of time series, for z of Eq. 0. L = 400. Window width 
W = 10. Triangles are simulational data. Thick line is fit 
to Eq. 0, with b = 0.924(2), oi = 1.03(1), a 2 = 0.51(1). 
X'd.o.f. = 1.5 x 10" 3 . 



these latter clearly differ from a pure Gaussian). 

We then attempted to fit the data in Fig. 0] to ana- 
lytical forms. The W — 10.0 results strongly suggest a 
double-Gaussian ansatz, as: 



*(«) = bG 1 (z) + (l-b) G 2 (z) 



(9) 



where Gi is a Gaussian centered at a* with variance of. 
As W grows, one would expect b — ► 1, a\ — ► 1 in Eq. J2j. 
Data for W = 10.0 are well fitted by 6 = 0.924(2), o x = 
1.03(1), a 2 = 0.51(1), as seen in Fig. El The \ 2 P er degree 



of freedom (x 



is 1.5 x 10 , indicating that the 



form Eq. lO indeed provides a satisfactory description of 
simulational results in this case. 

We have found that a similar fit, albeit of somewhat 



reduced quality (x 2 



3 x 10~ 3 , with b = 0.955(5), 



oi = 1.02(1), a 2 = 0.06(3)) is feasible for the W = 2.5 
data as well. Turning to W — 1.0, the double-Gaussian 
ansatz worked surprisingly well, producing x 2 d f = 6 x 
10~ 4 , with b = 0.53(5), a x = 1.24(4), a 2 = 0.77(1) (i.e. 
the two curves are roughly symmetric about z = 1, with 
approximately equal weights). 

Given that a double-peak structure is far from obvious 
for the W = 1.0 data, alternative forms must be con- 
sidered which might also provide a suitable fit to data 
in this limit. We investigated the family of roughness 
PDFs for l/f a noise 0, keeping in mind that win- 
dow boundary conditions (WT5C) are the appropriate 
ones in this case 0, 0, IE I^Hfill . Such PDFs are usually 
available in closed form [2||. However, close to a = 1 



Figure 7: (Color online) Crosses: scaled roughness distribu- 
tion ^(y) of time series, for y of Eq. (Jsj. L — 400, window 
width W = 1. Dashed line is double-Gaussian fit to data 
(Eq. 10). Full line is Fisher-Tippet-Gumbel distribution with 
window boundary conditions. Vertical axis is linear in (a), 
logarithmic in (b). 



it is more time-efficient to evaluate PDFs numerically 
via the usual procedure of first generating a very long 
sequence of Gaussian white noise, Fourier-transforming 
that sequence, multiplying the Fourier components by 
j-q/2 an£ j then inverting the Fourier transform 0, EI- 
The resulting sequence is pure 1// Q noise, which is then 
chopped into windows for analysis of the corresponding 
roughness PDF. 

The best fit of the l// a family to our data was 
achieved for a = 1, that is, the Fisher-Tippet-Gumbel 
(FTG) statistics of extremes 0. Even so, significant dis- 
crepancies remain. The overall picture is illustrated in 
Fig. where we have switched again to the variable y 
of Eq. (|HJl because the FTG curve is better visualized in 
this way 0,3- One sees that, even though the double- 
Gaussian curve gives an excellent fit in the central area 
of the plot where ^(y) > 0.1, it fails away from there, es- 
pecially at the lower end. As to the FTG curve, while it 
follows the data closely, it never actually matches them. 



V. DISCUSSION AND CONCLUSIONS 

The usual approach to the frequency domain in BN lit- 
erature is via the study of power spectra [Hi. l2^|. It has 
been found that, in the adiabatic limit of the inter- 
face model under consideration here, the power spectrum 
behaves approximately as l// 2 within an intermediate 
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range of frequencies. One might construe this as indi- 
cating that the pure l// 2 noise model of a Wiener pro- 
cess 0,0] applies in this case. However, the numerically- 
obtained full roughness PDF, which contains much more 
information than a section of the power spectrum, tells a 
more nuanced story. Indeed, in general it does not follow 
a shape close to that of 1/ f a curves, except for narrow 
windows. Even there, the closest fit within that family is 
for a ~ 1. 

The question then arises of whether the generalized 
Gaussian approximation underlying l// a noise models, 
in which the Fourier modes are considered as uncorre- 
lated |4], is suitable for the description of BN time se- 
ries. Our results, when considered in their evolution as 
window width varies, appear more consistent with the 
idea that the similarity of our PDFs to that of 1// noise, 
found at the narrow-window limit, is fortuitous. We re- 
call that, even in studies of real-space interface roughness, 
it is known that the independent-mode approximation 
gives rise to small but systematic discrepancies against 
experimental data, which can be traced back to higher 
cumulants of the correlation functions Q. Furthermore, 
even more severe discrepancies have been found when 
boundary conditions other than periodic (e.g. window, 
as is the case here) are considered pj.l27l|. 

Turning now to the double-Gaussian picture, admit- 
tedly phenomenological in its inspiration, nonetheless it 
gives a description which is both numerically closer to 
actual data, and spans a broad range of window widths. 

The physical origins of the double-peak structure may 
be traced back to the demagnetization term, and the con- 
sequent negative correlations illustrated in Fig.0 A win- 
dow of width W contains at least W/ tm segments whose 
internal roughness profiles are uncorrelated to each other. 
On the other hand, within each such segment, negative 
correlations are significant at least to some extent, thus 
preventing fluctuations from becoming very large. This 
latter effect gives rise to the secondary peaks at y < 0, 
or equivalently, z < 1. With tl ~ 0.14, tm — 0.5 for the 
L = 400 systems which have been the focus of our study, 
one has for W = 1 that both inter- and intra-segment 
fluctuations have similar weights, hence the b ~ 0.5 result 
for the double-Gaussian fit in that case. For W/tm 3> 1 
the dominant picture is one of many uncorrelated "blobs" 
of length ~ tm, yielding the effective single-Gaussian 
limit observed. 

The double-Gaussian picture displays features which 
are not fully understood at present. Fig. |H1 exhibits 
the variation of parameters b, a\, and a-i of Eq. 
against W for not very large window widths (in addi- 
tion to W = 1.0, 2.5, 5.0 and 10.0 we ran simulations 
at W = 0.5 and 1.5). While b varies approximately as 
expected within this theoretical framework (albeit with 
small non-monotonicities), and a\ follows a rather mono- 
tonic trend, the behavior of ai is intriguing, showing an 
apparent trend reversal. So far we have not able to pro- 
vide an explanation for this. 

An alternative explanation for the observed behavior 
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Figure 8: Fitting parameters 6 (triangles), ai (squares), and 
ai crosses) of double-Gaussian ansatz of Eq. against win- 
dow width W, for W = 0.5, 1.0, 1.5, 2.5, 5.0, 10. L = 400. 



at W ~ 1 may be proposed, following a line similar to 
that advanced for the evolution of r\ = data with in- 
creasing W (see Fig. and the respective discussion). 
In this scenario, the W ~ 1 PDF shapes would be non- 
universal (i.e. neither l/f a nor double-Gaussian). For 
larger W < 10 the central limit theorem would imply 
that, for the superposition of many (almost) decorre- 
lated non-universal profiles, effective Gaussian structures 
should emerge. In this view, the peak at z < 1 would 
again be ascribed to segments within which negative cor- 
relations are felt, with the peak at larger z corresponding 
to inter-segment profiles. 

Whatever the explanation of the behavior of rough- 
ness PDFs for W ~ 1, the extent of window widths for 
which an effectively double-peaked structure shows up is 
considerably larger than, say, T£, Thus, a fairly straight- 
forward way to detect the presence of demagnetization 
effects in experimental setups would be via the analy- 
sis of roughness PDFs of the induced signal V(t). Con- 
sidering e.g. the conditions for the Perminvar samples 
described in Ref. 0, where the average spacing between 
peaks is 13 msec and tm — 200 msec, analysis of windows 
of width ~ 2 sec should produce a well-defined double- 
peaked structure similar to that of Fig. H3 
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